On the critical behavior of disordered quantum magnets: The relevance of rare regions 
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The effects of quenched disorder on the critical properties of itinerant quantum antiferromagnets 
and ferromagnets are considered. Particular attention is paid to locally ordered spatial regions that 
are formed in the presence of quenched disorder even when the bulk system is still in the para- 
magnetic phase. These rare regions or local moments are reflected in the existence of spatially 
inhomogeneous saddle points of the Landau-Ginzburg- Wilson functional. We derive an effective 
theory that takes into account small fluctuations around all of these saddle points. The resulting 
free energy functional contains a new term in addition to those obtained within the conventional 
perturbative approach, and it comprises what would be considered non-perturbative effects within 
the latter. A renormalization group analysis shows that in the case of antiferromagnets, the previ- 
ously found critical fixed point is unstable with respect to this new term, and that no stable critical 
fixed point exists at one-loop order. This is contrasted with the case of itinerant ferromagnets, 
where we find that the previously found critical behavior is unaffected by the rare regions due to an 
effective long-ranged interaction between the order parameter fluctuations. 
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I. INTRODUCTION 

The influence of static or quenched disorder on the 
critical properties of a system near a continuous phase 
transition is a very interesting problem in statistical me- 
chanics. While it was initially suspected ±hat quenched 
disorder always destroys any critical point ,□ this. was soon 
found to not necessarily be the case. Harriso found a 
convenient criterion for the stability of a given critical 
behavior with respect to quenched disorder: If the corre- 
lation length exponent v obeys the inequality v > 2/D, 
with D the spatial dimensionality of the system, then 
the critical behavior is unaffected by the disorder. In the 
opposite case, v < 2/D, the disorder modifies the criti- 
cal behavior n This modification may either (i) lead to a 
new critical point that has a correlation length exponent 
v > 2/D and is thus stable, or (ii) lead to an uncon- 
ventional critical point where the usual classification in 
terms of power-law critical exponents looses its meaning, 
or (iii) lead to the destruction of a sharp phase transi- 
tion. The first possibility is realized in thej-conventional 
theory of random- T c classical ferromagnetsjj and the sec- 
ond one is probably realized in classical ferromagnets in 
a random field.cl The third one has occasionally been at- 
tributed to the exactly solved McCoy- Wu modelJa This 
is misleading, however, as has recently been emphasized 



in Ref. ^; there is a sharp, albeit unorthodox, transition 
in that model, and it thus belongs to category (ii). 

Independent of the question of if and how the criti- 
cal behavior is affected, disorder leads to very interesting 
phenomena as a phase transition is approached. Disorder 
in general decreases the critical temperature T c from its 
clean value Tj?. In the temperature region T c < T < T° 
the system does not display global order, but in an in- 
finite system one will find arbitrarily large regions that 
are devoid of impurities, and hence show local order, with 
a small but non-zero probability that usually decreases 
exponentially with the size of the region. These static 
disorder fluctuations are known as 'rare regions', and the 
order parameter fluctuations induced by them as 'local 
moments' or 'instantons'. Since they are weakly coupled, 
and flipping them requires to change the order parame- 
ter in a whole region, the local moments have very slow 
dynamics. GriffithsEI was the first to show that they lead 
to a non-analytic free energy everywhere in the region 
T c < T < Tj? s which is known as the Griffiths phase, or, 
more appropriately, the Griffiths region. In generic clas- 
sical systems this is a weak effect, since the singularity 
in the free energy is only an essential one. An impor- 
tant exception is the McCoy- Wu modelJa which is a 2-D 
Ising model with bonds that are random along one di- 
rection, but identical along the second direction. The 
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resulting infinite-range correlation of the disorder in one 
direction leads to very strong effects. As the temperature 
is lowered through the Griffiths region, the local moments 
cause the divergence of an increasing number of higher 
order susceptibilities, x (n) = d n M/8B n (n > 2), with M 
the order parameter and B the field conjugate to it, start- 
ing with large n. Even the average susceptibility proper, 
diverges at a temperature T x > T c , although the 
average order parameter does not become non-zero until 
the temperature reaches T c . This is caused by rare fluctu- 
ations in the susceptibility distribution, which dominate 
the average susceptibility and make it very different from 
the typical or most probable one. 

Surprisingly little is known about the influence of the 
Griffiths region and related phenomena on the critical 
behavior. Recent worka on a random- T c classical Ising 
model has suggested that it can be profound, even in 
this simple model where the conventional theory predicts 
standard power-law critical behavior, albeit with critical 
exponents that are different from the clean case. The au- 
thors of Ref. H have shown that the conventional theory 
is unstable with respect to perturbations that break the 
replica symmetry. By approximately taking into account 
the rare regions, which are neglected in the conventional 
theory, they found a new term in the action that actu- 
ally induces such perturbations. In some systems replica 
symmetry breaking is believed to be associated with ac- 
tivated, i.e. non-power law, critical behavior. Reference 
^ thus raised the interesting possibility that, as a result 
of rare-region effects, the random- T c classical Ising model 
shows activated critical behavior, as is believed to be the 
case for the random-field classical Ising model,cl although 
in the case of the random- T c model no final conclusion 
about the fate of the transition could be reached. 

Griffiths regions also occur in the case of quantum 
phase transitions (QPTs), i.e. transitions that occur at 
zero temperatv«e-as a function of some non-thermal con- 
trol parameter Jj't3 Their consequences for the critical be- 
havior are even less well investigated than in the classical 
case, with the remarkable exception of certain 1-D sys- 
tems. Fishertj has investigated quantum Ising spin chains 
in a transverse random field. These systems are closely 
related to the classical McCoy- Wu model, with time in 
the quantum case playing the role of the 'ordered direc- 
tion' in the latter. He has found activated critical be- 
havior due to rare regkuis. This has been confiriped by 
numerical simulationsilJ Other recent simulationsE3 sug- 
gest that this type of behavior may not be restricted to 
1-D systems, raising the possibility that exotic critical 
behavior dominated by rare regions may be generic in 
quenched disordered quantum systems, independent of 
the dimensionality and possibly also of the type of disor- 
der. 

In this paper we consider this problem analytically for 
two QPTs in D > 1. We first concentrate on a simple 
model for a quantum antiferromagnet. Previous work,E3 
which did not take into account rare regions, had found 
a transition with some surprising properties. One of our 



goals is to check whether these results survive taking into 
account rare regions. We find that they do not; the pre- 
viously found critical fixed point is unstable with respect 
to the rare regions, and one finds runaway flow in all 
of physically accessible parameter space. We will discuss 
possible interpretations of this result. We then show that 
the critical behavior of itinerant quantum /erromagnets 
is not affected by the rare regions, in sharp contrast to 
the antiferromagnetic case. A brief report of some of our 
results has been given previously in Ref. 

The paper is organized as follows. In Sec. || we derive 
an effective action for an itinerant antiferromagnet in the 
presence of rare regions. In Sec. Ill we perform a one-loop 



renormalization group analysis of this action, and show 
that there is no stable critical fixed point to that order. 
In Sec. IV we perform an analogous analysis for itinerant 
ferromagnets and show that the previously found critical 
fixed point is stable with respect to the rare region effects. 
In Sec. [v| we discuss our results. Various technical points 
are relegated to three appendices. 



II. AN EFFECTIVE ACTION FOR DISORDERED 
ANTIFERROMAGNETS 

A. The model 

Our starting point is Hertz's actionEl for an itiner- 
ant quantum antiferromagnet. It is a -theory with 
a p -component order parameter field (p whose expecta- 
tion value is proportional to the staggered magnetization. 
The bare two-point vertex function is 



r (q, u n ) =t + q 2 



(2.1) 



with to the mean distance from the mean-field critical 
point, q is the wavevector, and uj n denotes a bosonic 
Matsubara frequency. We measure q and uj n in suitable 
microscopic units to make them dimensionless. As in Ref. 
|l3| , we modify this action by adding disorder in the form 
of a 'random-mass' or 'random-temperature' term. That 
is, we add to to a random function of position, 5t(x), 
which obeys a distribution with zero mean and variance 
A, 



(at(x)) = o , 



(<Jt(x)at(y)>=A5(x-y) 



(2.2a) 



(2.2b) 



For the sake of simplicity, we have taken the distribution 
to be delta-correlated. The two-point vertex now reads 

T(x - y, r - t') = r (x - y, r - r') 

+5(x - y) S(t - t') St(x) . (2.3) 

Here r denotes imaginary time, and To (x, r) is the 
Fourier transform of Fo(q, ui n ) in Eq. fl2.ip . The action 
reads 
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S[(/ ) } = S G [(t>]+u / dx {4>{x) ■ <t>{x)) 



(2.4a) 



with the Gaussian part given by 

So \ < t>] = \J dx d V 4>(?) T{x-y) <P(y) . (2.4b) 

Here we have introduced a four-vector notation, x = 
(x, t), J dx = J dx Jq^ T dr, and we use units such that 
K = Icb = 1. 

At this point, the conventional procedure would be 
to integrate £>ut the quenched disorder by means of the 
replica trick.! This would lead to an effective action that 
does not contain the disorder explicitly any longer, and 
that therefore does not easily allow for saddle-point so- 
lutions that are not spatially homogeneous. While the 
effective action would still be exact, this latter property 
would make it hard to incorporate the physics we are 
concentrating on in this paper. We will therefore take 
a different approac h, an d consider saddle-point solutions 
of the model, Eqs. (2.4), before integrating out the disor- 
der. Qur procedure roughly follows the one by Dotsenko 
et alJ3 for classical magnets. As we will see, however, 
there are important differences between the classical and 
quantum cases. 



B. Saddle-point solutions 

Let us consider saddle-point solutions of Eqs. ( [2.4|) that 
are time independent. For simplicity, we also consider a 
scalar field, p = 1, that we denote by </>(x). It will be 
obvious how to generalize the following considerations to 
the case p > 1. With these restrictions, the saddle-point 
equation reads 

(t + <ft(x) - dl) sp (x) + 4u0 s 3 p (x) = . (2.5) 

Although 4p } (x) ee is of course always a solution, in- 
homogeneous solutions also exist provided that-5f (x) has 
'troughs' that are sufficiently wide and deep.liH In Ap- 
pendix |A] we demonstrate this for a one-dimensional toy 
problem. We have solved Eq. ( |2.5| ) numerically for ro- 
tationally invariant potentials <5i(x) = /(|x|), and have 
found behavior that is qualitatively the same as in the 
one-dimensional model. Thus, if St(x) has one sufficie ntly 
deep and wide trough, there will be a solution of Eq. ( |2.5| ) 
that is exponentially small everywhere except within the 
trough region, where it shows a single hump. There are 
actually two equivalent single- hump solutions, one pos- 
itive, and the other negative. We denote the positive 
one by ^^(x) ee ^W(x). As is intuitively obvious, and 
demonstrated in Appendix |a|, it leads to a lower free 
energy than the homogeneous solution 0sp*(x) = 0. Al- 
though in principle any saddle point can be used as the 
starting point for a loop expansion, it is therefore reason- 
able to assume that the inhomogeneous one will already 



incorporate physics that would be much harder to obtain 
if we expanded about the homogeneous saddle point. 

Next consider a potential <$i(x) that contains many 
troughs that support an essentially non-zero local order 
parameter field. This will result in a saddle-point solu- 
tion that contains many regions of local order, which we 
will refer to as 'islands'. Of course, for an arb itrary po- 
tential St(x) it is not possible to solve Eq. (2.5) in closed 
form. However, as long as the concentration of the is- 
lands is low, as will always be the case sufficiently deep 
in the disordered phase (i.e., for sufficiently large to), the 
values of tp sp outside of the islands will still be expo- 
nentially small. If St has troughs leading to TV islands, 
which individually would result in positive saddle-point 
solutions tpi + \x), (i — 1, . . . ,N), it is therefore a rea- 
sonable approximation to write the global saddle-point 
solution as a linear superposition of the V| > For inde- 
pendent islands, there are actually 2 N equivalent saddle 
points, which we write as 

N 



»)(x) = $W(x)=£<r?^+>(x) 

N 

= J>(x) , (2.6) 



where a = 1, ...,2^ numbers the equivalent saddle 
points, and the af are random numbers whose values 
are either +1 or — 1. They thus obey a probability dis- 
tribution 



with 



(2.7a) 



(2.7b) 



Alternatively, one can consider the ipi{~x) random func- 
tions that are equal to either plus or minus (x) . 

For later reference, let us briefly discuss the effects of 
a weak interaction between the islands as a result of the 
exponentially small overlap between the functions ^>i(x) 
centered on different islands. One effect will be that the 
total amount of the order parameter on each island will 
not necessarily be equal to plus or minus the amount re- 
sulting from that island being ideally ordered, but that 
small deviations from this amount will be possible. If 
we still assume that the islands are statistically indepen- 
dent, we can model this effect by using a pr obab ility dis- 
tribution for the <Ji that is given by Eq. (2.7a) with a 
distribution p(o~) that is a broadended versio n of t he bi- 
modal delta-distribut ion it (a) given in Eq. (2.7b). We 
thus generalize Eqs. ( |2 . 7a|j2 . 7b|) to 



p[{<}\=I[pW) 



(2.7c) 



3 



For our purposes we will not need to specify the distribu- 
tion p(a) explicitly. It will turn out that an interaction 
between the islands, no matter how small, leads to new 
physics compared to a model where these interactions are 
neglected. 

We also note that the islands will have some dynamics, 
both due to interactions between the islands and due to 
interactions between an island and its immediate neigh- 
borhood. In principle, one could try to build this effect 
into the saddle-point approximation by looking for time 
dependent saddle points. However, this dynamics is ex- 
pected to be very slow due to the inertia of the islands. 
Moreover, the zero frequency component in a frequency 
expansion is expected to yield the dominant effect. We 
therefore restrict ourselves to static saddle points. 



C. The partition function for a given disorder 
realization 

Of the 2 N saddle-point solutions (x) discussed in 
the previous subsection, let us pick one, say to ex- 
pand about: 

(j)(x) = $W(x) +tp(x) . (2.8a) 
Then the partition function can be written 

Z[6t(x)] = [ D[<p(x)]e- s l* W W + *W' st W , (2.8b) 



where we show <5i(x) explicitly as an argument to em- 
phasize that we are still working with a fixed realiza- 
tion of the disorder. Equation ( |2.8b ) is exact as long 
as the integral extends over all fluctuations <p(x) of the 
field configuration. However, in practice the integral over 
ip{x) cannot be performed exactly, and in a perturbative 
treatment one restricts oneself to small deviations ip{x) 
from the chosen saddle point. Typical pairs of saddle 
points picked from the 2 N $w represent field configura- 
tions that are globally very different. They will thus be 
far apart in configuration space, with large energy barri- 
ers be tween them. (We will justify this statement in Sec. 
V A 2 below.) Expanding about one of the saddle points, 
as in Eqs. (fj7g|), is therefore not expected to yield a good 
representation of the partition func tion if one evaluates 
the functional integral in Eq. ( 2.8b ) perturbatively. On 
the other hand, the same argument suggests that we can 
simply sum the contributions to Z obtained by expand- 
ing about all of the 2 N saddle points, provided that we 
restrict ourselves to small fluctuations about each saddle 
point, 



2 N 

Z[5t(x)] « £ / D[y{ 

a=l J < 



x e 



-S[4> (a) (x)+ V (:r),<5t(x)] 



(2.9) 



amounts to an arithmetic average over the perturbative 
contributions coming from the vicinities of all saddle 
points. This average is our approximative way of tak- 
ing into account non-perturbative effects. 

Since we are interested in the effects of fluctuations 
about the saddle points, we next subtract .the saddle 
point action from the exponent in Eq. (2.9).t3 That is, 
we write 

Z[St(x)]a^2 f D[<p(x)} e -AS[»W(x) >(P (»),«(x)] j (210) 

a J< 

where 

AS*[$( Q \ tp, St] = 5[$ (a) + if, St] - 5[$ (a) , St] 

= S[p{x)]+4u ( dxv 3 {x)<$> {a \x) 



(2.11) 



So far we have implicitly assumed that there is no in- 
teraction between the islands. In reality, there will be a 
small interaction, one effect of which wi ll be to replace 
the bimodal distribution, Eqs. ( [2.7a 2.7b ), by the broad- 
ended distri bution given in Eq. (2.7c). The sum over a 
in Eq. ( 2.1C ) is then replaced by an integral over the of, 
weighted by the distribution p(cr). The partition function 
can now be written as 



+6u / dxip 2 (x)<f> (a \x) 2 



Z[6t(x)} « / D[ V (x)\ 



-{S[ V (x)]+5S[ V {x)]) 



with the correction to the action, SS, given by 



(2.12a) 



-8S[ip(x)] 



N 



Here J indicates an integration over small fluctuations 
only. Apart from a normalization factor, this procedure 



i=l j 
x e -4u j dx<p 3 (x) 2.D-«V t C+) (x) 

x e -6u J dx<p a (x) £. jCTi <T 3 ^ +) (x)^ +) (x) 

(2.12b) 

As mentioned previously, it is crucial to incorporate a 
small interaction between the islands. Indeed, if we used 
the distribution function, Eqs. ( |2.7[ ), for the er,; that is 
appropriate for non-in teract ing islands, then we could do 
the a- integral in Eq. ( 2.12b ) exactly. As we will show in 
Sec. Ill, the resulting action would not lead to new phys- 
ical effects compared to Ref. [l^. Apart from the broad- 
ening of the distribution, there are other, similar effects 
of the island- island interaction that we will neglect. For 
instance, in the second exponent on the right-hand-side 
of Eq. ( |2.12b ), the absence of any overlap between ip^ 

and ipj^ for i ^ j makes the spatial integral in that term 
vanishes unless i — j. This is no longer true for inter- 
acting islands. However, we will neglect this effect and 
still ta ke this term to be proportional to 5ij. Equation 
(|2.12t ) can then be written 
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s e -6«/^cp 2 ( a; )(^+)(x))^ I , (212c) 



Here (. . .) denotes an average over the a with respect to 
the broadened bimodal distribution p(a). The er-average 
is carried out by means of a cumulant expansion in pow- 
ers of our order parameter fluctuations <p(x). To the or- 
der (p A we obtain 

6S[<p(x)] =6u f dx^^Y, (^ (+) (x)) 2 (o?) 

+Au f dx^ 3 (x)J24 + \ x ) (oi) 
— 18u 2 J dxdyip 2 (x) tp 2 {y) 

x^(^ +) (x)) 2 (v4 +) (y)) 2 ((^)-<^ 2 ) 2 ) • 

i 

(2.13) 

Now (cri) = 0, (erf) > and (of) - (ct 2 ) 2 eee c, > 0. The 
last relation is only valid for a broadened distribution 
p(cr) which arises from interactions between the islands. 
For the original distribution 7r(cr), c, = and so the 
0(</? 4 ) term would vanish. If we collect all contributions 
to the action for one particular disorder configuration, 
we obtain 

S[<p(x)] + 8S[<p(x)] 

= \l dxdy (p{x) T(x - y) y{y) +u j dx tp 4 (x) 
+ 18u 2 



dxdy p 2 (x) V 2 (y) £ Cl (x) V^ (+) (y )) • 

i 

(2.14) 



Here we have used the fact that the first term in ( 2.13| ) 
only renormalizes the random-mass term in the Gaus- 
sian action. We will show in Sec. Ill that truncating the 



action at 0(</? 4 ) is justified since all higher order terms 
are irrelevant (in a power-counting sense) with respect to 
both the Gaussian fixed point and the antifcrromagnetic 
fixed point found in Ref. n3[ 



D. The effective action 

So far we have considered one particular realization of 
the disorder. In order to derive an effective action we 
now need to perform the disorder average. It is impor- 
tant to remember that the Landau functional, Eq. (2.14), 
depends on the disorder in two different ways: explicitly 
through the random mass in the Gaussian action, and 
implicitly through the saddle-point solutions ipf (x) that 
depend on St(x). 



The quenched disorder average over 5t(x) is performed 
via the replica trickjH which is based on the identity 



{log Z} St = lim 



{Z n }st - 1 



n— »0 n 



(2.15) 



Here {. . .}st denotes the disorder average. This results 
in an effective action S e s which is defined by 

{z n } st =/n D[ V a {x)\ { e -£ Q (^)]+w(*)])}^ 

a— 1 

n 

Y[D[cp a (x)]e- s °«& a W . (2.16) 



In the absence of 5S, carrying out the disorder average 
yields the usual terms that are familiar from the conven- 
tional theory. Up to 0(ip 4 ) they are: 

t;Y1 J dxd y ^°( x ) T o{x-y)^ a (y) 

+ uY,f dx( V a (x)f 

- A W dxdrdT* ((^(x,r)) 2 (/(x,r')) 2 . (2.17) 

Taking into account the additional term, 6S[(p a (x)], is 
more subtle since the functions ipi( x ) are implicit func- 
tions of St(x). We handle this problem by means of a 
cumulant expansion. To lowest order, the contribution 
of SS to the effective action is just the disorder average 
of SS, 

{SS}st =w f dxdy cp 2 (x) <p 2 (y) D^fx.y) , (2.18a) 



where 



and the correlation function 



D 



(2) 



( x > y) = { c > 



(2.18b) 



Si 



essentially describes the probability for x and y to belong 
to the same island. The properties of these correlation 
functions depend on the precise nature of the disorder. If 
the microscopic disorder <$t(x) is short-range correlated, 
as we have assumed in our model, then the island size 
distribution will generically fall off exponentially for large 

l(2) 



sizes. 



In this case the correlation function D> sl ' is also 



short-ranged in space. Keeping only the leading term in 
a gradient expansion, we can then replace it by a spatial 
(5-function. The case of an island size distribution that 
has a power-law tail (e.g. due to long-range correlations 
in the microscopic disorder) is discussed in Appendix 

Collecting all contributions to the effective action S e g 
up to 0(ip A ), absorbing a constant into w, and restoring 
the vector nature of the order parameter field, we finally 
obtain 
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Seff[v a (z)] = ~E / dxd y T (x~y) ( p a (x)-<p a (y) 
+U Y1 f d*dT{<p a (x,T) ■ <p Q (x,r)) 2 
— /J (A + w 5 a p) / dxdrdr' 

x (<p«(x, r) • <^(x, r)) (</(x, r') • </(x, r')) . 

(2.19) 

The ui-term is generated by taking into account the in- 
homogeneous saddle points. A perturbative expansion 
about the homogeneous saddle point, as was performed 
in Ref. [l], misses this term. It has the time structure 
of the random-mass or A-term, and the replica structure 
of the quantum fluctuation or u-term. We will further 



discuss its physical meaning in Sec. V A. In the following 



section we will show that the critical fixed point found in 
Rcf. Ul is unstable with respect to this new term in the 
action. 



III. RENORMALIZATION GROUP ANALYSIS 

A. Tree-level analysis 

Let us firstiustify our truncation of the Landau expan- 
sion in Sec. [lj by showing that all terms of higher than 
quartic order in ip are irrelevant (in the renormalization 
group sense) by power counting with respect to the crit- 
ical fixed point of Ref. 13 To this end, wc analyze the 
effective action, S e g, Eq. (2.1i), at tree level. 



Let us denote the scale dimension of any quantity Q 
by [Q] , and define the scale dimension of a length L to be 
[L] = — 1. The scale dimension of the imaginary time is 
[t] = — z, which defines the dynamical critical exponent 
z. We first analyze the Gaussian fixed point. From the 
s tru cture of the two-point vertex function Tq given in Eq. 
(2.1), we see that u> n scales like q 2 . This implies z = 2. 
The scale dimension of the field <p is found (from the 
requirement that the action must be dimensionless) to be 
[y(a;)] = D/2. The scale di mensio ns of the coefficients of 
the terms of 0(<p A ) in Eq. ( 2.19 ) are found to be [u] — 
2 — D, and [A] = [w] = 4 — D. Thus, u is irrelevant with 
respect to the Gaussian fixed point as long as D > 2, 
while A and w are relevant for D < 4. The Gaussian fixed 
point is therefore unstable, and we will have to perform 
a loop expansion close to D = 4 in the next subsection. 

We now show that all terms of 0(tp 6 ) and higher are 
irrelevant with respect to the Gaussian fixed point. First 
of all, there are the conventional terms of the schematic 
form 



U2m / dxp 2m (x) 



(3.1) 



with coupling constants u 2m (v,i = u). These are irrel- 
evant since [u 2m ] = 2— (m — l)D < 0. In addition to 



these terms, the cumulant expansion of ( [2. 12c ) generates 
higher order terms with more time integrations than the 
conventional terms for a given power of ip. For instance, 
at 0(ip 6 ) we have two terms, 



8u 2 J dx dy ^ (x) v 3 (y)J2 ^ +) « ^ +) W ^ 



(3.2a) 



and 



-36u 3 y dxdydz <p 2 (x) ip 2 (y)ip 2 (z)^2(jp ( i + \^)) 

i 

x(V, (+) (y)) 2 (V 4 (+) (z)) 2 

x((^)-3(a.f)<0+6<0 3 ) . (3.2b) 
Upon averaging over the disorder these terms become 

v e J dxi dx 2 ^{xi) ip 3 (x 2 ) C^(xi,x 2 ) , (3.3a) 

and 



(3.3b) 



w 6 I dxidx 2 dx 3 ip 2 (xi)ip 2 (x 2 )(p 2 (x 3 ) 

xD isl ( X X,X 2) X3) , 



respectively, with vq oc u 2 , wq oc u 3 . The correlation 
functions c£j (x, y) and dS^(x, y, z) are defined analo- 
gously to D$ (x, y) in Eq. fl2.18b|) , and are related to the 
probability for x, y and x, y, z, respectively, to belong to 
the same island. We again concentrate on the generic 
case where the island size distribution falls of exponen- 
tially for large islands (for a discussion of other cases, 
see Appendix |b|) . In this case both correlation functions 
are short-ranged and can be localized for power-counting 
purpose s. Th is eff ective ly leaves only one spatial integral 
in Eqs. (B.3a) and ( 3.3b| ). Therefore, the scale dimensions 
of the coefficients are [ve] = 2(2 — D) and [wq] = 2(3 — D). 
Consequently, both terms are irrelevant with respect to 
the Gaussian fixed point near D = 4. 



More generally, we obtain from Eq. (2.12c) terms that 
contain powers of y 3 , terms that contain powers of ip 2 , 
and mixed terms that contain both ip 3 and ip 2 . For power 
counting purposes, the most relevant term for a fixed 
power of ip is the one with the most time integrations. 
For even powers of ip, these are the terms 

w 2m J dxi . . .dx m p 2 (x 1 ) . . .p 2 {x m ) D^(xi, . . . ,x m ). 

(3.4) 

Localizing the correlation function _D( m ) we find for 
the scale dimension of the coupling constant [w2m] = 
2m — (m — 1)D. Terms with odd powers of <p are always 
less relevant than the preceding term of even order. We 
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FIG. 1. The three diagram struc ture s that contribute to 
the flow equations shown in Eqs. (3.5). The dashed lines 
stand for any of the vertices whose coupling constants are u, 
A, or w. 



conclude that all terms of higher than quartic order are 
irrelevant with respect to the Gaussian fixed point near 
D = 4. 

So far we have determined the scale dimensions with 
respect to the Gaussian fixed point. At the non-trivial 
critical fixed point discussed in Ref. O the anomalous 
dimension of the field ip is r\ = + O(t) since, as in the 
ordinary theory, there is no wavefunction renormaliza- 
tion at 1-loop order. This implies that all results on the 
irrelevancy of the terms of order ip e and higher carry over 
from the Gaussian fixed point to the non-trivial critical 
fixed point found in Ref. O. 



B. Perturbation theory, and flow equations 

In the last subsection we have shown that the Gaussian 
fixed point is unstable for D < 4. We must therefore 
carry out a loop expansion for the effective action, Eq. 
19). To control the perturbation theory we considp: 
4 — e spatial dimension and e T time dimensions. t£l 



(1 



D 

This leads to the replacement of / rfr by J drr^" 1 . At 
the Gaussian fixed point the scale dimension of the field 
tp is now [tp] — (d + ze T — 2)/2. In the same vein the scale 
dimensions of it, A, and w are [it] = e — ze T , [A] = e, and 
[w] = e, respectively. The perturbation theory becomes 
a double expansion in e and e T . 

To obtain the renormalization group flow equations, w% 
perform a frequency-momentum shell RG procedure.LS 
The diagrams that contribute to the renormalization of 
the coupling constants u, A, and w are shown in Fig. [I]. 
To one-loop order, we obtain the following flow equations, 



^ = 0-2e r )u-4(> + 8W 2 +48uA , (3.5a) 
dl 



4- = eA + 32 A 2 - 8 (p + 2) u A + (2 A - w) w , 
dl T £t 



dw 
~dT 



ip 



< ir -h - w 2 - 8 (p + 2) u w + 24 (2 A -w)w 
1 £t 

+ 8w 2 . (3.5c) 



The mass t of the two-point vertex, which describes the 
distance from the critical point, is of course also renor- 
malized. However, since we are interested in the stability 
of a critical fixed point, it suffices to consider the flow 
on the critical surface. The factors of T~ €t in Eqs. ( |3.5| ) 
arise from the fact that some diagrams that contain the 
u>-vertex lead to Matsubara frequency sums without an 
accompanying temperature factor. Since the critical sur- 
face for the quantum phase transition is defined by T = 
in addition to t = 0, the natural coupling constant for the 
T = flow is w = wT~ e T. Putting T = 0^3 the flow 
equations can then be rewritten in the form 



du 



2e 7 



4 (p + 8) u 2 + 48 u A 



(3.6a) 



dA 
If 



dw 
~dl 



eA + 32A 2 -8(p + 2)uA + 8pAw; , (3.6b) 



(e - 2e T ) w + 4p w 2 - 8 (p + 2) u w + 48 A w 



(3.6c) 



C. Fixed points and their stability 



(3.5b) 



The flow equations, Eqs. ( |3.6|) , possess eight fixed 
points. The fixed-point values of the coupling constants, 
and the corresponding eigenvalues of the linearized RG 
transformation are listed in Table |. Four of the fixed 
points (Nos. 1-4 in Table Q) have a zero fixed-point value 
of w, w* — 0. These are the fixed points studied before 
in Ref. [l3| The other four fixed points have w* ^ 0. 

Let us first consider fixed point No. 4. This is the 
critical fixed poin.t that was found within the conven- 
tional approach.EJ We find that the local moments, rep- 
resented by the u;-term, render this fixed point unstable 
for p < 4, since in this case the third eigenvalue, A^,, 
is positive. However, for p > 4 the ui-term is irrelevant 
with respect to this fixed point, and the fixed point is 
stable fox.4 < p < p c . To one-loop order, and for e = e T , 

Pc = iee 

A stability analysis for the new fixed points shows that 
they are all unstable for p < 4 with the exception of No. 
8. At this fixed point, w has a negative value, while Eq. 
( [2. 13 ), with reasonable assumptions about the distribu- 
tion p(a), yields a positive value for the bare value of w. 
Since the structure of the flow equations does not allow 
w to change sign, we conclude that this fixed point is 
unphysical for generic realizations of the disorder. It is 
interesting to note, however, that fixed point No. 8 is 
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FP values 
A* 



eigenvalues 

A2 At, 



e-2e T 
4(p+8) 



(e±er) 
16 (p-1) 






-e/32 

(4-p)e+4(p+2)e T 
64 (p-1) 



6 - 2e T 
-(e-2e T ) 
2e T - e/2 



8(p+4) £x -(p-4)6 
P+8 



F 3 ! 



P-1 



e - 2e T 

(p-4)(e-2e T ) 
p+8 

2e T - e/2 

(4-p) £ +4(4-p)e T ) 
4 (p-1) 





e — 2e T 
4 (P+8) 



c+4e T 
8(10-p) 






(4e T - e)/64 

((p-4) £ +24e T ) 
32 (10 -p) 



-(e-2e T )/4p 

(p-4)(<r-2<r,-) 
4p (p+8) 

-(4e T + e)/16p 



(p-4)( £ +46 T ) 
8p(10-p) 



e/4 - e T /2 
-(e-2e T ) 
e T + e/4 



4e T — e 

2(p+4)2e T -(p-4) £ 
P+8 



-c+Vc 2 -o 

p-1 



2(l()-p) 



-(e-2e T ) 

(p-4) (e-2 Q 
P+8 

2e T - e/2 

(p-4)(e+4e T ) 
2(10-p) 



TABLE I. Fixed points of the flow equations, Eqs. (3.6) and the eigenvalues of the corresponding linearized RG trans- 
formation, p is the number of order parameter components. Ai,2 are the eigenvalues in the absence of w, and X w is the 
additional eigenvalue. A, B, C, and D are defined as A = (3e — 4e T )p + 16e r , B = 16(p — 1) (e + 4e r ) [(4 — p)e + 4(p + 2)e T ], 
C= (16-p)e + 4(p-4)e T , and D= (10 - p) (p + 8) (e + 4e T ) [(p - 4)e + 24)e T ]. 



stable against replica symmetry breaking (see Appendix 
0)- 

For p < 4, and to one-loop order in our double expan- 
sion in powers of e and e T , there is thus no stable fixed 
point. Consistent wit h th is, a numerical solution of the 
flow equations, Eqs. (3.6), shows runaway flow in all of 
physical parameter space. We will discuss the physical 
meaning of this result in Sec. [v| below. 



IV. THE CASE OF ITINERANT 
FERROMAGNETS 

In Ref. |l5| a generalized LGW functional for the fer- 
romagnetic transition in a disordered itinerant electron 
system was derived starting from a fermionic description. 
The effects of rare regions were not explicitly considered 
in this work. Here we show that although the rare regions 
were neglected in the explicity calculations in that paper, 
the effective field theory derived in Ref. [l5] still contains 
these effects. We will further show that taking them into 
account does not change the previous conclusions. 

We first briefly recall the effective action that was 
derived in Ref. 15l In the long-wavelength and low- 
frequency limit, the replicated action is given by 

S«b,i = \^2 fd Xl dx 2 T (x 1 -x 2 )M a (x 1 )-M a (x 2 ) 



x(M a ( Xl ) ■ M a (x 2 )) (M a (x 3 ) ■ M a (x 4 )) 
" A 5l/ d * dTdT ' (M Q (x,r)) 2 (M^r')) 2 . 

(4.1) 



a,0 ' 



Here M is the order parameter field whose expectation 
value is the magnetization, and u and A are coupling con- 
stants that in general are wavenumber and frequency de- 
pendent. An important point is that these coupling con- 
stants in general do not exist in the limit of zero frequen- 
cies and wavenumbers, i.e. the effective action describes 
a non-local field theory. This is because in the process 
of deriving a LGW functional that depends only on the 
order parameter field, soft (viz., diffusive) fermionic de- 
grees of freedom have been integrated out. In writing Eq. 
(4.1) we have used that the coupling constant A is finite 
in the long- wavelength limit, so that it can be treated as 
a number. U4, on the other hand, is singular in this limit, 
see Eq. ( |4.3| ) below. For small wavenumbers the Fourier 
transform of the two-point vertex T$ is given by 



r (q, uj n ) = t + u 2 (q) 



with, 



Here u 



(u-2) 



,(D-2) 



and u?' are finite numbers. Note that for 



1 13-2 



U 



(°)„2 



(4.2a) 



(4.2b) 



D < 4 (in particular, in th e physical dimension D = 3), 
the first term in Eq. (4.2b) dominates the second one as 
q — > 0. U4, in wavenumber space at zero frequency, is 
schematically given by 
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u 4 (q -> o) 



iqr e +o(| q r 4 ) 



(4.3) 



i.e., 114 diverges for D < 6. The singularities in the 
wavenumber dependences of u 2 and M4 mean that the 
field theory is non-local. As mentioned above, their phys- 
ical origin are diffusive fermionic particle- hole excitations 
that were integrated out in deriving Eq. (4.1). 

Next we argue that in at least one well defined physical 
situation, it is easy to u nco ver the effects of rare regions 
that are implicit in Eq. ( 4T ) . The basic argument is that 
on length scales small compared to the elastic mean free 
path £, the field theory is effectively local. This implies 
that as long as the local moments or instantons decay on 
a scale A < £, the techniques discussed in Sec. |l| can be 
used to include the effects of these inhomogeneous saddle 
points on the final long-wavelength theory. Further we 
will show that the quenched randomness that leads to 
these local moments is implicitly contained already in 
the last term in Eq. (4.1). To this end we first note, cf. 
Ref. [l^ and below, that the exact critical behavior near 
t he f erromagnetic transition can be determined from Eq. 
Second, we partially undo the replica trick in Eq. 
by writing the logarithm of the partition function, 



(hi 



or the free energy, as Eq. ( 2.15| ) with Z on the right hand 



side given by, 



Z = J D[M(x)] exp(-S[M(x),St(x)}) , (4.4a) 
with, 

S[M(x),8t(x)} = i J dxt dx 2 T(x t - x 2 ) 

xM(xi) • M(x 2 ) 

+ J dxi . . . dx^ u^ixx, X2, x$, X4) 

x(M(x 1 )-M(x 2 ))(M(x 3 )-M(x 4 )) . 

(4.4b) 



T(x) is given by 



r(x) = r Q (x) + st(x) 



(4.5) 



with To from Eq. (4.2a). <5t(x) in Eqs. (4.4) is a random 
function of position that is Gaussi an d istributed with the 
first two mome nt s giv en by Eqs. (|2.2[). It trivaliy follows 
that Eqs . ( 2.15 ), (4.3) - ( [4.5| ), and ( |2.2| ) are equivalent to 
Eq. Q. 

The arguments given in Sec. || for the AFM case apply 
equally well to the action given by Eq. (4.4b). To make 
this more precise let us consider the non-localities. In 
terms of scaled variables, Eq. (4.2b) for u 2 (q) can be 
written to lowest order in the disorder as, 



/ \ (0) 
M 2 (q) = u\ 



_q_ 

hp 



_c_ /|q| 



D-2 



(4.6) 



with kF the Fermi wavenumber, I the elastic mean free 
path, and c an interaction dependent constant that is at 



most of ord er unity. In general the expansion implied 
by Eq. ([0]) assumes that |q| <C kp. In terms of length 
scales let A ~ l/|q| be the scale over which the order 
parameter varies, and Xp ~ kp 1 the Fermi wavelength. 
Further, to be specific we consider the physical dimen sion 
D = 3. The analytic, square gradient, term in Eq. (4.6) 
then dominates the seco nd t erm when A -C £■ That is, 
the non-locality in Eq. (4.6) is irrelevant when spatial 
scales shorter than a mean free path are considered. For 
Xp <C A <§; I we then have, 



"2(q) 



,(0) 



— y 

kp) 



(4.7) 



Note that for this argument to be valid we need Xf /£ ~ 
X/kpi <C 1, i.e ., weak disorder is required. Similarly, 
(i4 in Eq. (4.3) can be replaced by a constant when 
Xf <C A -C I. The net result is that when the local 
moments vary on a length scale smaller than £, they can 
be described by a local field theory analogous to the one 
discussed in Sec. || even though the long-wavelength the- 
ory is non-local. If we assume, as we did in the antifer- 
romagnetic case in the previous section, that the island 
distribution falls off exponentially for large island sizes, 
this will always be true for sufficiently small disorder. 

With the above ideas and the techniques developed in 
Sec. [n| the final long wavelength theory to descibe the 
ferromagnetic phase transition, explicitly including the 
effects of rare regions, is determined by the action, 



Scff = S c s,i + SS, 



off: 



(4.8a) 



with i§ e ff,i given by Eq. (4.1) and 



SS cff = -wJ2 f dxdrdr 1 (M a (x, t) • M a (x, r)) 

x(M Q (x,r') -M q (x,t')) , (4.8b) 

where w is a finite constant. 

Power counting immediately reveals that the coupling 
constant w, just like A, is an irrelevant operator with 
respect to the Gaussian fixed point discussed in Ref. [l^. 
The rare regions therefore do not change the critical be- 
havior in this case. The physical reason for this is the 
effective long-range interaction between the order param- 
eter fluctu ations that is described by the |q| D ~ 2 -term in 
Eq. (4.2b). This stabilizes the Gaussian fixed point by 



suppressing all fluctuations, including the static disorder 
fluctuations resonsible for the local moments. 



V. DISCUSSION AND CONCLUSION 

In this section, we conclude by discussing the results 
obtained in the previous sections. 







A. General considerations 

We begin our discussion by considering the physical 
underpinnings of some general aspects of our technical 
procedure. 



1. Local moments and annealed disorder 

Let us first of all give a simple physical interpreta- 
tion of the u>-term in the effective action, Eq. ( [2.19 ) 



which is the most important of the contributions that 
reflect the existence of rare regions and local moments. 
Since the local moments are self-generated by the elec- 
tronic system, in response to the potential created by the 
quenched disorder, they are an integral part of the sys- 
tem and in equilibrium with all other degrees of freedom. 
In our approximation, which takes into account only the 
static local moment fluctuations, the effect of the rare 
regions therefore amounts to the existence of static, an- 
nealed disord er. Indeed, a straightforward generalization 
of Eq. ( 2.10 ) is to integrate over a manifold of saddle 
points 3?(x), weighted with an appropriate distribution 
P[*(x)], 



J £>[*(x)]P[*(x)] J D[ V (x)] e 



(5.1) 



which makes obvious the annealed-disorder character of 
the average over the saddle points. The detailed result 
of the integration over the saddle points will of course 
depend on the distribution P, which in turn depends on 
the microscopic details of the disorder realization in the 
system. However, any physically reasonable distribution 
will lead in particular to a term in the effe ctive action 
that has the structure of the ui-term in Eq. ( [2.19 ). Since 
the saddle points are separated by large energy barri- 
ers in configuration space (see Sec. VA2j below), this 
term clearly cannot be obtained by perturbatively ex- 
panding about the trivial homogeneous saddle point as 
is done in the conventional theory. Thus, our method 
approximately takes into account what one would call 
'non-perturbative' effects in the usual approach. 

It is im porta nt to note that the new term in the ac- 
tion, Eq. ( 2.19 ), differs from the usual quantum fluctu- 
ation or w-term only in its time structure. In the clas- 
sical limit, therefore, w just renormalizes u, decreasing 
its bare value. This is indeed well known to be the only 
effect of static annealed disorder in classical systems. In 
their analysis of classical magnets, the authors of Ref. || 
therefore considered a more elabor ate s cheme for doing 
the su m ov er saddle points in Eq. ( 2.1C ), or the integral 
in Eq. (5.1), that leads to a term that breaks the replica 
symmetry. Our way of approximating that integral can 
be considered as a zeroth order step in the approxima- 
tion scheme of Ref. @. In the quantum case, the time 



structure results in this zeroth step already giving a non- 
trivial result, and in this sense quantum systems are more 
sensitive to rare region effects than classical ones. The 
physical meaning of replica symmetry breaking in this 
context is not quite clear. However, in the quantum case 
it is not necessary to enter into this discussion. The AFM 
fixed point is unstable already under the effects consid- 
ered above, and no new fixed point exists. Considering 
replica symmetry breaking in addition to our effect would 
not change this conclusion. In the FM case, it turns out 
that the previously found Gaussian fixed point is stable 
against replica symmetry breaking as well as against the 
quantum effect, as we will discuss in more detail below. 



2. Energy barriers between saddle points 



A question that arises in connection with Eq. (2.9) or 
(|5.l|) is whether it is really true that there are large energy 
barriers between the various saddle-point configurations, 
as our approximation for the partition crucially depends 
on this assumption. Let us first consider the case of the 
Ising model (p = 1), for which we performed the explicit 
derivation in Sec. O. Suppose we have two saddle points 
that differ only by the sign of the order parameter on 
one particular island. In order to turn one of these spin 
configurations into the other, we need to flip all of the 
spins on that island. (For simplicity, we refer to the order 
parameter field as 'spins'.) In order to do so, one must go 
through an intermediate state with a domain wall across 
the island. The energy of that domain wall can be esti- 
mated from the squared gradient term in the free energy, 
integrated over the island, J J dx (V0(x)) 2 , with J the 
coupling between the spins. The thickness of the domain 
wall is a microscopic length a, and hence the energy of the 
domain wall, or the energy barrier between the two sad- 
dle points, is proportional to J L D ^ 1 a/a 2 = J L D_1 /a, 
with L the linear size of the island. In the case of a 
continuous spin model (p > 1) an analogous argument 
holds, except that now all length scales are of order LE 
This leads to JL D " 2 for the energy of a domain wall. 
For D > 2, there is thus only a quantitative difference 
between the Ising case and the continuous spin case. 

In either model, the domain wall energy will have to be 
multiplied by the number of islands by which two typical 
saddle points differ. For the Ising case, let us consider N 
islands, with 2 W saddle points and 2 N ~ 1 {2 N — 1) pairs of 
saddle points. The probability distribution {p_/v(n)}, for 
a pair of saddle points to have n islands that are different 
is easily found to be 



p N {n) 



1 



2 N -1 



(5.2a) 



For large N, this becomes a Gaussian distribution with 
mean N/2 and variance y/N /2, 



PN^eo(n) = 



V2^V 



,-2(n-N/2) 2 /N 



(5.2b) 
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One expects this to be true for the continuous spin case 
as well, although the statistical analysis becomes much 
more involved in that case. The miscroscopic energy of a 
domain wall thus gets multiplied by a macroscopic num- 
ber, leading to energy barriers between almost all pairs 
of saddle points that go to infinity in the thermodynamic 
limit. This justifies our approximation. 

Finally, we note that our considerations maximize, and 
probably overestimate, the effects of local moments or 
disorder induced instantons. The discussion above seems 
to imply an exponential number of saddle point solutions 
that are unrelated by symmetries, with barriers between 
them that approach infinity in the bulk limit. This in 
turn implies an exponential number of thermodynamic 
states, or a finite complexity. Such a proposition is con- 
troversial in other contexts, e.g. for spin glasses. How- 
ever, in our considerations we have effectively neglected 
the interactions between the local moments. One antici- 
pates these interaction to correlate and weaken the rare 
regions. Indeed, in Ref. [23] it was argued that long-range 
interactions that arise from the itinerant nature of the 
electrons quench most of the local moments. If this hap- 
pens in the systems we consider, then we likely overesti- 
mate the number of distinct thermodynamic states. It is 
also possible that our theory is valid only in an interme- 
diate time region, and that the interactions between the 
local moments must be taken into account in the limit of 
asymptotically long times. 



3. Nature of the local-moment phase 

Another point we have not yet addressed is the physical 
nature of the phase that is induced by the presence of 
the local moments. In order to show that we are dealing 
with a Griffiths phase, let us consider the local moment 
contribution to the order parameter susceptibility, xlm- 
Let us adopt a ferromagnetic language for simplicity, and 
denote the magnetic moment on the island number i by 
AL. Then we have 



Xlm 



1 



l/T 



(Mi) (Mj)) 



(5.3) 



si. 



where (. . .) denotes a thermodynamic average. Since 
there is no overall magnetization, Y], (Mi) = 0, and in 
our saddle-point approximation the island magnetization 
is static. This yields 



Xlm 



e 



si 



const 
T 



(5.4) 



where the constant is given by {J2i(Mf) / J2i Vi}st- We 
see that the order parameter susceptibility diverges for 



T — > whenever there are islands, and in our sim- 
ple saddle-point approximation the divergence takes the 
form of a Curie law. Our saddle point thus really de- 
scribes a Griffiths phase. 



4- Finiteness of the free energy 

As a final gen eral p oint, let us again consider the effec- 
tive action, Eq. (2.19), which determines the free energy. 
Since the w-term and the it-term have the same struc- 
ture except for an extra time integral in the former, it 
seems as if the w-term contributes a term to the free en- 
ergy that diverges as the temperature g oes to zero. One 
has to keep in mind, however, that Eq. (2. IE) represents 
a Landau expansion that has been truncated at 0(ip 4 ). 
It is easy to see that higher order terms in the Landau 
expansion lead to even more str ongl y divergent contribu- 
tions to the free energy, see Eq. (3.4). This simply means 
that the loop expansion for the free energy of a quantum 
system with static annealed disorder is singular, and a 
resummation to all orders would be necessary to obtain 
a finite result. From a RG point of view, which holds 
that the higher order terms in the Landau expansion are 
irrelevant, the solution of this paradox lies in the fact 
that, if a fixed point existed, it would be w that has a 
finite fixed-point value, not w. Since w — wT (for the 
physical case e T = 1), this ensures that the fixed-point 
Hamiltonian has a finite free energy. 



B. Results for the AFM case 



As we have shown in Sec. Ill, and reiterated above 



taking into account the rare regions in the AFM case de- 
stroys the stability of the fixed point found in Ref. |l3| , 
and one finds runaway flow in all of the physically acces- 
sible parameter space. Three possible interpretations of 
this result are, (i) there is no transition to a state with 
long-range order, (ii) there is a transition, but the corre- 
sponding fixed point is inaccessible by perturbative RG 
techniques, or (iii) there is a fluctuation-induced first or- 
der transition (which causes the runaway flow). The last 
conjecture can be checked by calculating the free energy 
to one-loop order and then explicitly verifying whether it 
has a double minimum structure as a function of the q&. 
der parameter. We have performed such a calculation.Eil 
and found that this is not the case. This rules out sce- 
nario (iii). 

On the basis of our results, we cannot decide between 
scenarios (i) and (ii). Scenario (i) would imply that arbi- 
trarily weak disorder necessarily destroys quantum AFM 
long-range order. This is an unlikely proposition, but it 
cannot be ruled out at present. The alternative is sce- 
nario (ii), i.e. the existence of a non-perturbative fixed 
point. The nature of such a fixed point, if it exists, is 
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a priori unclear. The analogies with l-D systems men- 
tioned in the Introduction, as well as Ref. |l2|, suggest 
that an unconventional infinite disorder fixed point with 
activated scaling is a possible interpretation of the run- 
away flow. However, there also could be a conventional 
fixed point that is not accessible by our methods. In 
this context it is intere sting to note that the case p > 4 
discussed in Sec. Ill C provides an example of a stable 
conventional fixed point that describes a transition with 
power-law critical behavior in the presence of rare re- 
gions. 

Let us also come back to the fact that in Sec 
found a stable fixed point (No. 8 in Table 



[II 



we 
?0. 



with w 

As was pointed out in Sec. Ill C , for generic realizations of 
the disorder, which lead to a positive bare value of w, this 
fixed point is unphysical since it has w* < 0. However, 
mathematically one can have w < for certain choices of 
the distri bution P[{erf }] in Sec. || that are more general 
than Eq. ( 2.7c ). This leaves open the possibility that at 
least in some systems there is a stable, conventional crit- 
ical fixed point that is accessible with our method. We 
note that this fixed point is stable against replica sym- 
metry breaking, see Appendix This is in contrast to 
the case of classical magnets,El where all fixed points are 
unstable against replica symmetry .breaking, and remi- 
niscent of the result of Read et alE3 on quantum spin 
glasses, where the quantum model was also found to be 
more stable against replica symmetry breaking than its 
classical counterpart. The technical reason for this en- 
hanced stability i s ver y similar to the point discussed 
at the end of Sec. VA, namely that the parameter that 
would induce replica symmetry breaking appears as T 
times a finite fixed point value, and hence vanishes at 
the quantum critical point. 



C. Results for the FM case 

For itinerant quantum ferromagnets, we have found, 
that the rare regions do not affect our previous results. E3 
The physical reason for this is the long-range interactions 
between the spin fluctuations in these systems. They are 
induced by soft modes in the itinerant electron system 
and stabilize the Gaussian critical behavior against fluc- 
tuations, including the static disorder fluctuations that 
lead to local-moment formation. A crucial point for our 
conclusion is the survival of these long-range interactions 
in the presence of local moments, so it is worth discussing 
this in some detail. I— . 

The derivation of the long-range interaction^ shows 
that its origin is soft spin-triplet particle-hole excitations 
in the electron system. An obvious question is whether 
local moments act effectively as magnetic impurities that 
give these soft modes a mass. If this were the case, 
then the singular wavenumber depe nden ces |q| d_2 and 
|w„|/q 2 in the Gaussian vertex, Eq. ( 4.2a ), would be cut 
off and the ferromagnetic effective action would have the 



same structure as the antiferromagnetic one. The answer 
to this question is not obvious since the local moments 
are self-generated, and thinking about them as analo- 
gous to externally introduced magnetic moments can be 
misleading. This is underscored by the fact that the rare- 
regions/local-moment physics enters the theory in the 
form of annealed disorder, as we have seen in Sec. VA 



above. In our derivation of the effective action, Sec. IV, 
the wavenumber singularities are not cut off. The physics 
behind this is that both singularities are consequences 
of spin diffusion, which in turn is a consequence of the 
spin conservation law. The rare regions ultimately derive 
from a spin-independent disorder potential, which clearly 
cannot destroy spin conservation. The long-range inter- 
actions between the spin fluctutations are therefore still 
present in the bare effective action, and hence the Gaus- 
sian fixed point is stable in our tree-level analysis. We 
note, however, that at present we cannot rule out the pos- 
sibility that loop corrections might lead to qualitatively 
new terms in the action. If such new terms included a 
RG-generated spin dependent potential, then this might 
change our conclusions. This would not necessarily vi- 
olate the spin conservation argument given above, since 
an effective spin-dependent potential acting only on the 
itinerant electrons, which are not taking part in the lo- 
cal moment formation, could change the critical behav- 
ior while preserving spin conservation for the system as 
a whole. 

A more detailed investigation of this point is not fea- 
sible within the existing framework of the ferromagnetic 
theoryEJ This is because in the existing theory all degrees 
of freedom other than the order parameter, including var- 
ious soft modes, have been integrated out. This leads to 
a non-local field theory which is unsuitable for an explicit 
loop expansion. A remedy would be to derive an effec- 
tive theory that keeps all soft modes explicitly and treats 
them on equal footing, leading to well-behaved vertices 
that allow for explicit calculations. This project is left 
for future work. We also note that at such a level of 
the analysis one should also include effects due to inter- 
actions between the rare regions, which we have mostly 
neglected. Such interactions are known to weaken the 
effects of the rare region,E3 but in general it is not known 
by how much. 



D. Summary, and Outlook 

In summary, we have studied the effects of rare dis- 
order fluctuations, and the resulting local moments, on 
itinerant ferromagnets and antiferromagnets. Techni- 
cally, this has been achieved by considering non-trivial 
saddle-point solutions before performing the disorder av- 
erage. A perturbative RG analysis of the resulting effec- 
tive field theory incorporates effects that would require 
non-perturbative methods within a more standard pro- 
cedure. In the ferromagnetic case we have found that, 
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at least within our level of-analysis, the previously found 
quantum critical behaviorLI is stable with respect to local 
moment physics. The reason is an effective long-range 
interaction between the spin fluctuations that strongly 
suppresses fluctuations, stabilizing a Gaussian critical 
fixed point. In the antiferromagnetic case, however, we 
have found that the local mopents destroy the previ- 
ously found critical fixed pointO To one-loop order, and 
for order parameter dimensionalities less than 4, no new 
fixed point exists and one finds runaway flow in all of 
physical parameter space. This may indicate either the 
absence of long-range order, or a transition that is not 
perturbatively accessible within our theory. 

An important technical conclusion is that for quantum 
phase transitions, and within the framework of a repli- 
cated theory, rare regions can have a qualitative effect 
already at the level of a replica-symmetric theory, in con- 
trast to the case of classical magnets.B The ferromagnetic 
fixed point, which was found to be stable against the 
replica-symmetric quantum effects induced by the rare 
regions, is also stable against replica-symmetry breaking. 

We have concentrated on the role of fluctuations about 
a non-trivial, but fairly crudely constructed, saddle-point 
solution of the held theory. It would also be interesting to 
study a somewhat more sophisticated saddle-point theory 
in more detail, and to determine the detailed properties 
of the Griffiths phase in such an approximation. 

Finally, we mention that our methods are not specific 
to magnets, and can be applied to other quantum phase 
transitions as well. For instance, it is believed that for a 
complete understanding of the properties of doped semi- 
conductors, and of the metal- insulator transitions ob- 
served in such systems, Jt_is_necessary to consider the 
effects of local moments. oEKj This can be studied with 
the methods developed in this paper. 
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W e fir st consider the one-dimensional counterpart of 
Eq. 

(t + St(x)-d 2 x ) <f>{x)+g4?{x)=Q , (Al) 

with a simple square well potential, 

St M = f-V for0<x<a (A2) 
w 10 elsewhere v ' 



Standard methods lead to a solution inside the well (0 < 

x < a) 



/zJT cni/i^ x 2 

(/>in[x) = J— - V , 

V 9 dn^/v 2 x 2 



(A3a) 



where cn an dn are elliptic functions, and a solution out- 
side of the well, 



to 2\[2ce 



Here 



with 



^2,1 



± vAx 2 -4/3 



a = 2(t - Vb) 



(A3b) 



(A4a) 



(A4b) 



c and [3 are constants of integration that are determined 
by the requirement that the solution and its derivative 
be continuous at x = a. Other solutions exist, but the 
one given is the only one that satisfies physical bound- 
ary conditions. Furthermore, the physical solution exists 
only for 



< [3 < a 2 /4 



(A5) 



To demonstrate the existence of a threshold, we expand 
the above solution for small values of v^a- To leading 
order in this small parameter, we obtain for the constants 
of integration 



v\ e % 



(A6a) 



and 



APPENDIX A: A ONE-DIMENSIONAL 
SADDLE-POINT EQUATION 

In this Appendix we discuss the saddle-point equation, 
Eq. (2.5), for a particular realization of the disorder po- 
tential <5i(x). In particular, we aim to show that the 
existence of a non-zero solution requires the width and 
depth of the potential well to be above a threshold, and 
that the non-zero solution lowers the free energy with 
respect to the zero one. 



R = _ _ J°_ 

' A 2a 2 



(A6b) 



From the condition for the existence of the physical so- 
lution, Eq. (|A5|), we see that a 2 a 2 > 2tQ is a necessary 
and sufficient condition for the solution to exist, which is 
the desired threshold property. 

The free energy in saddle-point approximation is sim- 
ply given by the saddle-point action. It is physically 
plausible that the non-homogeneous solution constructed 
above leads to a negative free energy, and is thus energet- 
ically favorable compared to the zero solution. We have 
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FIG. 2. t(r) (dotted line) and <j>(r) (solid line) for a = 1, b = 1, c = 1. 



ascertained this numerically for a large variety of well pa- 
rameter values, and have found the inhomogeneous solu- 
tion to lead to a negative free energy whenever it exists. 

T o sol ve the three-dimensional saddle-point equation, 
Eq. (|2.5|) , is much harder. For sperically symmetric wells, 
to + St(x) — t(r), with r = |x|, analytic solutions can still 
be found in closed form for special forms of the potential. 
By scaling <j> and t, the equation can be written 



V 2 + t(r)<p(r) +(f> 3 {r) = 



(A7) 



It is easy to show that for 
-26 



t(r) 



46< 



2 c 



(1 + cr 2 ) (1 + cr 2 ) 



(1 + br) 

4 b cr 8 c 2 r 2 

(l + cr 2 )(l + 6r) + (1 + cr 2 ) 2 
4c a 2 (l + 6r) 2 e- 2fcr 



2 b 2 



(1 



(l + 6r 2 ) 



(1 + cr 2 ) 2 



the physical solution is given by 

a(l + br) exp(— br) 



(r) 



(1 + cr 2 ) 



br) 
(A8a) 

(A8b) 



Here a, 6, and c are parameters that determine the shape 
of the well. In this case, the physical solution exists for all 
real values of the three parameters, but the form of the 
potential is such that the volume of the well cannot be 
smaller than some minimum value. This is the three- 
dimensional analog of the threshold behavior demon- 
strated above for the one-dimensional case. We have also 



solved the ODE, Eq. (A7), numerically for more general 
potential wells, and have found the same type of thresh- 
old behavior. 

As in the 1-D case, physical arguments suggest, and 
numerical integration confirms, that the inhomogeneous 



solution leads to a lower free energy than the homoge- 
neous one whenever the former exists. In Fig. || we show 
t he solution and the corresponding potential well, Eqs. 
( |AS| ), as a representative example of a locally ordered 
region in a 3-D system. 



APPENDIX B: ISLAND SIZE DISTRIBUTIONS 
WITH A POWER-LAW TAIL 



In Sees. Ill and IV we have assumed that the island 



correlation functions and C^T' are short-range cor- 

related, i.e. have a scale dimension of (m — 1)1?. Here 
we briefly discuss the extent to which one can relax this 
condition without changing our results. 

Suppose that the island-size distribution is power-law 
correlated, leading to scale dimensions of the above cor- 
relation functions that are given by (m — 1)(D — a) with 
a > 0. Let us consider the FM case first. The least 
irrelevant term, viz. Eq. ( gjtj ), remains irrelevant with 
respect to the Gaussian fixed point of Ref. [l5| as long as 
a < 4 - D (for 2 < D < 4). The D-dependence of this 
result reflects the fact that for D > 4 the effective interac- 
tion ceases to be long-ranged, and an ordinary mean-field 
fixed point is stable. All higher order terms in the action 
are less relevant than the ui-term. 

In the AFM case, the w-term is relevant with respect 
to the conventional fixed point even for a = 0. By power 
counting, we find the condition that none of the higher 
order terms become relevant as well, viz. a < D — 3 for 
D close to 4. Here the D-dependence reflects the fact 
that the coupling constant wq is marginal in D = 3 even 
for a = 0, see Sec. [II A. 
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APPENDIX C: STABILITY UNDER REPLICA 
SYMMETRY BREAKING 

In this Appendix we briefly consider the effects of 
replica symmetr y bre aking (RSB). A generalization of 
our action, Eq. ( 2.19 ), analogous to Ref. || that allows 
for RSB is 

Srt[<p a (x)] = ±Y, [ dxdy t p a (x)-r (x-y)<p a (y) 

a J 

+u]T f dxrfr(^(x,r)-( ¥3 Q (x,r)) 2 
— ^^(A + w Q/ 3) / dxdrdr' 

a,/3 J 

x(^(x,r) • V Q (x,r)) (</W) • (</(x,r')) . 

(CI) 

In Sec. [fi] we had w a p = w, which resulted in a replica 
symmetric theory. Now we. allow for 1-step RSB in 
Parisi's hierarchical scheme^ where w a p in the replica 
limit is parameterized by means of a step function with 
a parameter xq, 



w(x) = | 



w for < x < 1 
Wi for Xq < x < 1 



(C2) 



Defining w = wT Ct as before, and w = wi,T £t , we ob- 
tain the 1-loop flow equations 

^ = (e-2e r )u-4(p + 8)u 2 + 48uA , (C3a) 
dl 



dA 



e A + 32 A 2 - 8 (p + 2) u A + 8p A w 
8pxo Aw + 8p Aw , 



(C3b) 



dw 
~dT 



(e - 2 e T ) w + 4p w 2 - 8 (p + 2) u w + 48 A w 



-4p(l-x )w 2 , (C3c) 



(e - 2 e r ) w + 48 A u) - 8 (p + 2) u w 
+ &pww + 8p{l-x )w 2 . (C3d) 



The replica symmetric case is recovered by putting xq = 
1. We now perform a linear stability analysis of fixed 
point No. 8, with fixed-point values of u, A and w as 
given in Table |, and w* = 0. The first three eigenvalues 
are as shown in Table |, and the fourth one is A^, = 
X w = (p - 4)(e + 4e T )/2(10 - p). Fixed point No. 8 is 
therefore stable against 1-step RSB. Although the fixed 
point is unphysical for generic realizations of the disorder, 
as discussed in Sec, [II, this is an interesting contrast to 
the classical case,El where all fixed points are unstable 
against successive terms in the hierarchical RSB scheme. 
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